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ABSTRACT 

We develop a simple, time-dependent Comptonization model to probe the 
origins of spectral variability in accreting neutron star systems. In the model, 
soft "seed photons" are injected into a corona of hot electrons, where they are 
Compton upscattered before escaping as hard X-rays. The model describes how 
the hard X-ray spectrum varies when the properties of either the soft photon 
source or the Comptonizing medium undergo small oscillations. Observations of 
the resulting spectral modulations can determine whether the variability is due 
to (i) oscillations in the injection of seed photons, (ii) oscillations in the coronal 
electron density, or (iii) oscillations in the coronal energy dissipation rate. 
Identifying the origin of spectral variability should help clarify how the corona 
operates and its relation to the accretion disk. It will also help in finding the 
mechanisms underlying the various quasi-periodic oscillations (QPO) observed 
in the X-ray outputs of many accreting neutron star and black hole systems. As 
a sample application of our model, we analyze a kilohertz QPO observed in the 
atoll source 4U 1608-52. We find that the QPO is driven predominantly by an 
oscillation in the electron density of the Comptonizing gas. 

Subject headings: accretion, accretion disks — radiation mechanisms: thermal 
— stars: neutron — X-rays: stars 
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1. Accretion Disk Coronae 

The X-ray spectra of accreting neutron stars and black holes in binary systems 
typically appear to have at least two components (for recent reviews, see White, Nagase, & 
Parmer 1995; and Tanaka & Lewin 1995). The soft spectral component resembles blackbody 
radiation. It is thought to come from optically thick gas in the accretion disk and possibly, 
in the case of an accreting neutron star, from gas near the stellar surface. The hard spectral 
component comes from hotter gas that fails to reach thermal equilibrium with its radiation; 
the radiation intensity falls off with increasing photon energy roughly as a power-law, 
Ie oc up to a cutoff photon energy Ec beyond which the intensity drops exponentially. 
In the most rapidly accreting neutron star systems, the Z-sources, the "hard component" 
cuts off at energies ^lOkeV and appears at least partially thermalized; the hard component 
extends to higher energies in the atoll sources, which have lower accretion rates, and in 
black hole candidate systems. The multicomponent appearance of their X-ray spectra 
suggests that the accretion flows around these objects may be divided into a minimum of 
two substructures, one consisting of cooler, dense gas that for convenience we shall call 
the "disk" component, and a coronal component of hotter, tenuous gas. Observations 
of rapid variability in X-ray binary systems (see van der Klis 1995 for a review) also 
support this picture of discrete physical substructures in the accretion flow. Fluctuations 
in the hard spectral component tend to be much more pronounced than those in the soft 
component, and to be more weakly correlated with fluctuations in the soft component on 
short timescales than on longer timescales. 

Comptonization has long been recognised as the probable source of the hard spectral 
component in cosmic X-ray sources (Thorne & Price 1975; Shapiro, Lightman, & Eardley 
1976), since it naturally produces high-energy spectra of approximately the right shape, 
and is the dominant interaction between photons and electrons in hot, tenuous gases. (Here 
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we do not consider neutron stars with surface magnetic fields Bg ^ 10^^ Gauss, as their 
extremely strong magnetic fields make the physics of their X-ray spectra unlike that of 
other accreting neutron stars.) The hot gas cools by Compton scattering low-energy "seed 
photons" to higher energies, and the Compton upscattered photons escape the system as 
the hard spectral component. The source of the seed photons is unclear at present, and is 
likely to be different in black hole and neutron star systems. In some models for spectral 
formation in black hole systems (see, for example, Haardt & Maraschi 1993), soft photons 
from the optically thick disk seed the Comptonization, but in others (e.g., Narayan, Yi, 
& Mahadevan 1995) the soft photons are due to cyclotron/synchrotron emission from 
the hot corona itself. In neutron star systems, cyclotron emission from the neutron star 
magnetosphere can provide the seed photons, and it is probable that Compton emission 
from hot magnetospheric gas dominates the hard spectrum (Psaltis, Lamb, & Miller 1995); 
nevertheless, the disk, magnetospheric boundary layer, and stellar surface also may play 
important roles. The question of where seed photons come from depends on the geometrical 
relation between the hot and cold fiow substructures, which is poorly known at best. For 
example, in the case of accretion by a black hole, advection-dominated disk models (see 
Narayan 1996 and references therein) and the classic Shapiro, Lightman, and Eardley 
scenario have a cold outer disk which becomes hot and optically thin to absorption within 
some critical radius, so that the corona occupies the most central regions of the fiow, while 
the cold disk lies beyond it; however, models in which the coronal gas sandwiches the 
colder disk have also received considerable theoretical development (Liang & Price 1977; 
Bisnovatyi-Kogan & Blinnikov 1977; and Galeev, Rosner, & Vaiana 1979; more recent 
studies emphasizing spectral formation include Haardt & Maraschi 1993; Haardt, Maraschi, 
& Ghisellini 1997; and Dove, Wilms, & Begelman 1997). Given the uncertainties, even very 
general information about the coronal structure and soft photon source would be welcome, 
in black hole candidate and in neutron star systems. 
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Rapid X-ray variability is now thought to hold vital clues to the nature of the accretion 
flow "engines" in these systems. In particular, quasi-periodic oscillations have attracted 
attention from both observers and theorists, since they are believed to contain otherwise 
inaccessible information about the structure and operation of the central accretion flow. (A 
quasi-periodic oscillation, or QPO, is enhanced variability in a relatively narrow range of 
frequencies around a so-called "centroid" or QPO frequency, which thus resembles but is not 
identical to periodic modulation in the emission from a source.) The QPO so far discovered 
in X-ray binaries span a wide range of frequencies, from i^gpo ^ 1 Hz to over a kilohertz, 
and appear to belong to a number of phenomenologically distinct categories, each of which 
presumably has a different mechanism. Proposed mechanisms for QPO tend to be situated 
within or near the radiation-generating central regions of their accretion flows, and include 
modulations of the flow at the inner edge of the disk or at the magnetospheric boundary 
(for the classic "magnetospheric beat frequency" models, see Alpar & Shaham 1985, and 
Lamb et al. 1985; also see Miller, Lamb, & Psaltis 1997 for a beat frequency model applied 
to the recently discovered kilohertz QPO), wave-hke oscillations in the inner disk (e.g., 
Chen & Taam 1995, and references therein), and oscillations in a coronal flow component 
(Fortner, Lamb, & Miller 1989; Miller & Lamb 1992). Although leading candidate models 
have emerged for one or two of the QPO categories, satisfactory explanations exist for very 
httle of the rich QPO phenomenology that has come to light over the past decade. QPO 
remain promising, and enigmatic. 

Without assuming a detailed hydrodynamic model for the QPO, several attempts 
have been made to extract general information about the physical characteristics of the 
accretion flow and the QPO mechanism from the spectral structure of the observed QPO; 
i.e., how the amplitude of the QPO varies from one band of photon energies to the next, 
and the degree to which the oscillation in a given band lags behind those in other bands. 
Early EXOSAT-based studies of the 20-50 Hz QPO in Z-sources (Hasinger 1987, van der 
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Klis 1987) showed that X-ray photons at higher energies lag behind those at lower energies 
by several milliseconds, an effect which immediately suggested the progression of photons 
from lower to higher energies in a Comptonizing corona. The discovery stimulated several 
groups to perform calculations of time-dependent Comptonization in which the soft photon 
injection rate underwent small oscillations about its average value (Wijers, van Paradijs, & 
Lewin 1987; Stollman et al. 1987). In addition to injection rate oscillations, Stollman et 
al. (1987) examined the effects of small oscillations in the coronal electron density. They 
concluded that either type of modulation, in the soft photon source or in the coronal 
electron density, could account for the observed hard lag; the data available at the time 
were insufficient to discriminate between the two possibilities. 

Observations of accreting neutron star and black hole binary systems with the Rossi 
X-ray Timing Explorer (RXTE) satellite have revealed unprecedented detail in their X-ray 
variations, and motivate us to return to the issue of how much can be learned about 
X-ray coronae and the origins of X-ray variability through time-dependent Comptonization 
models. In this paper we develop a method of calculating spectral variations based on 
the time-dependent Kompaneets equation that is far more computationally efficient than 
the Monte Carlo methods employed by Wijers et al. and Stollman et al. Our calculations 
permit us to introduce small oscillations into any of the underlying properties of the 
Comptonization process and study how they drive the emerging radiation. Since the 
Kompaneets approach is best for systems with scattering optical depths r well above unity, 
our method is more appropriate for accreting neutron star systems (in which coronal 
optical depths are expected to be >5) than for black hole candidate systems (in which it 
is estimated r 1). Using values for coronal parameters expected to be characteristic of 
neutron star systems, we explore how X-ray spectral variations depend on coronal properties 
and on the nature of the underlying driving modulation. We show that considerable 
differences exist between the spectral variations produced by modulations in the soft photon 
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injection, in the coronal electron density, and in the coronal electron temperature. These 
differences should make it possible to fit observed spectral variations as a superposition of 
the different types of modulation, and to obtain constraints on QPO mechanisms and on 
system parameters such as the size of the Comptonizing corona that are more reliable than 
those obtained by considering only one type of modulation in isolation. 



2. Time-dependent Comptonization Model 

Since the properties of an accretion disk corona and its soft photon source are so poorly 
known, we adopt the simplest possible description for them consistent with our goal of 
physically interpreting variability in hard X-ray spectra. We ignore the effects of gradients 
in the electron density and temperature of the corona, and treat the Comptonization 
process as if it occurred in a completely homogeneous medium with electron density Uf. and 
temperature T. We do, however, allow these quantities to undergo small oscillations as 
photons in the corona are Comptonized. We further enforce the assumption that the corona 
is homogeneous by stipulating that once injected into the corona, any photon may escape 
from it with a probability per unit time 

Tesc = CaTUe/N^sc = l/(A^esc^c) , (1) 

where c is the speed of light and (Jt is the Thomson cross section. The escape rate Fesc 
depends only on the mean time between collisions, tc, and on the average number A^esc of 
collisions suffered by a photon before its escape. Thus the model omits any information 
about the location of a photon within the corona, and all of the properties of the corona are 
subsumed in the three quantities, T, tc, and A^esc- 

Our model also treats the injected photons simply. Photons enter the Comptonizing 
environment at a rate N. At the instant of its entry, the probability that a photon has an 



-8- 



energy between E and E-\-5E is given by fini{E)5E. In this paper, we will take the injected 
photon distribution to have a blackbody spectrum: fin]{E) = ATi^j^i?^(exp[£'/Tinj] — 1)^^, 
where ^4 is a normalisation constant. (We write the temperature in units of energy, so that 
Boltzmann's constant /cb = 1.) To study how fluctuations in the injection process influence 
the spectrum of photons escaping from the corona, we introduce small oscillations in N 
about its average value. We do not permit /inj to vary. Thus, we allow for variability in 
the overall rate of photon injection, but ignore more subtle modulations in the injected 
spectrum. 

For nonrelativistic coronal temperatures (T <^ mgC^) and sufficiently low photon 
energies {E <^ mgC^), the Kompaneets equation (Kompaneets 1957; see Psaltis and Lamb 
1997 for a complete derivation of the equation and a discussion of its limitations), 

tAf = -\dE [-^TEf + E-'f + TdE {E^f)] + Winj - ^ , (2) 

describes the temporal development of the distribution / of photons in the Comptonizing 
medium. The Kompaneets equation is essentially a photon continuity equation in energy. 
Its natural boundary conditions are that the term in square brackets (the photon current) 
vanish as £^ ^ and as £^ — > oo. 

The spectrum of photons escaping from the Comptonizing medium is given by 

fout{E) = Tesc/ , (3) 

SO that modulations in the observed spectrum can be driven both by changes in the 
spectrum / of photons within the cloud and by changes in the rate Fesc at which photons 
escape from the cloud. 

To find how variations in the Comptonization process affect the escaping spectrum, we 
assume that a given parameter in our model executes small oscillations around its average 
value, and solve for the resulting behavior of the escaping spectrum to flrst order in the 
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oscillating perturbation. For example, to see how /out reacts to variations in the photon 
injection rate N, we would write the injection rate as the sum of a constant value Nq and a 
small oscillatory part: 

N = No + Nie-''^^ . (4) 
Then to first order in the perturbation, the escaping photon spectrum is given by 



fout{E) ~ /outo + e /( 



out 1 ) 



where /outo is the escaping spectrum in the limit of a vanishing perturbation. 



/out 1 — /l / [^c^esc] 



' 



(5) 
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and the perturbation amplitude /i of the coronal spectrum is given by the linearised 
Kompaneets equation, 

d 



\Nesc / mgC^ dE 



ATEh - E'h - T^E^h 



(7) 



In addition to perturbations in the soft photon input, we also consider two types of 
perturbation in the corona itself. Perturbations in the rate of energy dissipation in the 
corona drive variations in its temperature: 



T = To + Tie- 



-iujt 



The escaping photon spectrum is still given by eq. (|^), but /i is now determined by 



iuU ) /i = ^ 
mpC"^ dE 



-4E (Ti/o + To/i) + EVi + 4^E^ (Ti/o + Toh] 

dE 



(9) 



Perturbations in the supply of matter to the corona can cause its size / or electron 
density Ue to change, altering tc and A^^esc We assume that A^^esc is connected to the electron 
density rig through A^^esc = t(t + 1), where the scattering optical depth of the corona is 
defined to be r = arnj. Thus, 
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(10) 
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(This choice is motivated by simplicity; the exact form of the relationship between A^esc a-nd 
tc or Ue depends on geometrical assumptions, and on how the length scale / is defined. As 
an example, for escape from the center of a homogeneous Comptonizing sphere, one finds 
A^esc = 37r~^T[T + (4/3)] when I is defined as the radius R of the sphere, but if one instead 
defines / = SV^-^i? ^ 0.6R, the result is a relationship between A^esc and tc close to the 
one we assume.) Now the variation in the escaping photon spectrum is given by 

/outl = -/outO f §^ + + /l/ [«sc]o , (11) 

V-'Ve.sr.O InO/ 



where 



\Nesco ^ rUeC^ dE 



4TEh - E^h - T^E^h 



Nh^M + h^ (12) 



governs the perturbation in the coronal spectrum. Although in general we could allow the 
parameters and iYggc to vary independently, in this paper we will assume that the length 
scale I of the corona remains fixed (we choose I — 10^ cm as characteristic of circumstellar 
regions in neutron star systems), while the electron density undergoes small oscillations. 

This simple model for dynamic Comptonization allows us to study the different spectral 
variabilities that arise from modulations in the supply of coolant (the soft photons), the 
supply of heat, and the supply of matter to the corona. In the following section we will 
examine the characteristics of each type of variability, paying special attention to traits 
that might serve to distinguish one from another observationally. 



3. X-ray Variability: Numerical Results 

Each type of perturbation — in photon injection, in coronal heating, and in the electron 

density of the corona — produces a different characteristic variation in the spectrum of the 
escaping radiation. By identifying the superposition of perturbations that best matches the 
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observed spectral variations of a given source, we can learn more about the physics of X-ray 
coronae, and in particular about their sources of material and energy. In this section we 
solve numerically for the observable spectral variations produced by each of the basic types 
of perturbation under a range of conditions. We find that, as functions of photon energy, 
the rms relative amplitude 



of each type of spectral variation differ considerably from the other types, and so should be 
readily distinguishable. 

We illustrate the main differences between the three types of spectral variability with 
a typical model corona. Coronal electrons with an average temperature of Tq = 10 keV 
Comptonize 0.1 keV seed photons. The photons suffer on average A^esco = 100 collisions 
before escaping from the corona. All the perturbations have frequencies v = uo / (27r) = 40 Hz 
and relative amplitudes of lO'^; i.e., Ni/Nq = 10"^ Ti/Tq = 10"^ and A^esci/A^esco = 10"^. 

Figure 1 shows the relative amplitudes (bold lines) and phases 0{E) (lighter lines) 
produced by each perturbation. Oscillations in the seed photon injection rate cause 
spectral variations with an approximately energy-independent relative amplitude (bold short 
dashes); since the residence time of photons in the corona, tesc = 1/resc, is short compared 
to the oscillation period, the spectral response is quasi-static. To a good approximation, 
the emergent spectrum simply rises and falls with the photon injection rate. The phase 
lag (light short dashes) is small and has a minimum near the seed injection energy Tinj; it 
rises monotonically with energy above Tinj. The phase lag corresponds roughly to the time 
delay ~ tc{mf.c^ /T) required by photons to upscatter from low energies ~ Tjnj to higher 
energies ~ Tq (as may be estimated from the Kompaneets equation ; also see Pozdnyakov, 
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(13) 



and the phase 




(14) 
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Sobol, & Sunyaev 1983). Temperature-driven spectral variations (solid lines) are very 
different. They have a minimum at photon energies near STq and rise sharply thereafter. 
The associated phase variations also contrast strongly with those due to oscillations in 
the photon injection rate, and reflect the redistribution of photons from energies below 
~ 3Tq to higher energies that occurs when the coronal temperature increases. Oscillations 
in the corona's optical depth (long dashes) alternately overpopulate and depopulate the 
photon spectrum above Tq, resulting in a relative amplitude that rises with photon energy 
near and above Tq, but less steeply than one driven by coronal temperature oscillations. 
The spectrum below ~ Tq falls as the optical depth rises. Oscillations at energies above 
To are nearly out of phase with those below, except for a phase lag that is caused by the 
upscattering time delay. Thus, each type of oscillation has its own distinctive signature. 

In the following subsections, we systematically examine how each of the spectral 
oscillation signatures depends on the frequency of the oscillation, and on the average 
properties of the Comptonizing medium and soft photon injection process. The perturbed 
spectra depend most strongly on the mean coronal electron temperature Tq, and on the 
ratio of the photon residence time to the oscillation period, vtesc- We choose Tq to be 5, 10, 
or 15 keV, characteristic of coronae in accreting neutron star systems. Injection energies Tlnj 
are drawn from the set {0.1 keV, 0.5 keV, 1 keV}. Oscillation frequencies oi u = 40, 100, 
and 1000 Hz are used to generate the perturbed spectra. Escaping photons experience an 
average of A^esco = 10, 100, and 400 collisions in our calculations, so that the combination 
utesc ranges from ~ lO^^-lO'^, and the spectral response goes from the quasi-static limit to 
fully dynamic behavior. 
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3.1. Perturbations in the Photon Injection Rate 

Suppose that the soft photon source oscillates shghtly, driving small oscillations in the 
emergent spectrum. If the oscillation frequency v is small, z/tcsc ^ li the response of the 
spectrum is close to the quasi-static [y —>■ 0) limit. Thus, at all but the highest frequencies, 
injection modulations simply cause the brightness of the source to undergo small oscillations 
about its average value, with little change in spectral shape. The relative amplitude is 
therefore fiat, and the phase shifts are small (see Figs. 2a and 2b). Spectral oscillations at 
higher i/ have larger phase lags and are more strongly attenuated at energies away from Tlnj . 

The dependence of the perturbed spectrum on the coronal temperature Tq is illustrated 
in Fig. 2c. The relative amplitude varies little with changes in Tq. Because the diffusion of 
photons in energy progresses more rapidly in a hotter corona, the phase lags of the spectral 
variations are smaller for a hotter corona than for a cooler one. Figure 2d shows that 
the amplitude of the spectral oscillations above 1 keV has almost no dependence on the 
injection temperature. The phase lags are somewhat more sensitive to Tinj, although the 
dependence is still not strong (approximately logarithmic in To/Tinj). The phase delay is 
minimised at energies near Tjnj, since photons escaping at energies near Tjnj typically have 
spent less time in the corona than photons at more distant energies. 

3.2. Perturbations in Electron Temperature 

Oscillations in the coronal electron temperature force photons from lower to higher 
energies. The spectral oscillations are thus anticorrelated with the temperature oscillations 
at low energies, and correlated at higher energies (Fig. 3a). The energy at which the 
behavior changes from anticorrelation to correlation corresponds to a minimum in the 
relative amplitude. 
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When the Comptonized spectrum is very nonthermal (i.e., when the "y-parameter" 
y = ANcscoTo/{mc'^) -C 1; see Pozdnyakov et al. 1983), the minimum occurs at energies 
somewhat above Tinj, but considerably lower than Tq. Fig. 3b shows that the minimum 
moves toward an asymptotic position at STq in the case of more thermalized spectra, 
for which the Wien peak at 2Tq is more prominent (these correspond to y » 1). 
This follows from the fact that in the neighborhood of the Wien peak, we expect the 
time-averaged spectrum to be /o ~ CnoTm.E^TQ^e~^/^°, where Cnorm is a normalisation 
constant independent of Tq. The To-dependence follows from the shape of the thermal 
peak, and the need to maintain a constant photon number (constant normalisation) at 
different electron temperatures. In the quasi-static limit the relative spectral variation is 
then Tie-''^^d\nfo/dTo, and so the relative amplitude is ^ |(-3/To + E/T^)Ti\. The 
approach to the quasi-static, thermalized limit is seen in Fig. 3c. The minimum in the 
relative amplitude occurs at 6kcV ~ 1.2To for the relatively unthermalizcd case Tq = 5keV, 
and has reached 33keV w 2.2Tq in the more thermalized case Tq — 15keV. 

As Fig. 3d shows, the influence of the photon injection temperature on the spectral 
oscillations is minor. It is slightly easier for photons introduced at higher energies to form 
a Wien peak, all other things being equal. Consequently, cases with higher Tlnj are a little 
better thermalized than cases with lower Tinj , and have amplitude minima at slightly higher 
energies. 

3.3. Perturbations in Electron Density 

When the coronal electron density oscillates, the spectral variations at energies near 
Tinj are nearly anticorrelated with the driving oscillations. The anticorrelation occurs 
because an increase in optical depth both suppresses the rate at which photons escape from 
the cloud and allows photons to evolve to energies more distant from Tlnj. The latter effect 
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is dominant at energies E '^Tq, where the spectral oscillation is nearly in phase with the 
electron density oscillation, except for a small phase lag due mainly to the time required for 
photons to diffuse out of the corona. 

The dependence of the spectral variations on driving frequency is shown in Fig. 4a. 
Unsurprisingly, the relative amplitudes of the variations only begin to betray their frequency 
dependence when the oscillations are fastest, and the spectral response emerges from the 
quasi-static limit. Since the time lag between the responses in different energy channels 
has little dependence on the frequency of the driving oscillations in the quasi-static limit, 
the phase lag at a high photon energies grows approximately linearly with frequency. As 
the mean electron density increases, the time lag increases, along with the escape time tesc- 
This is seen in the phase lags of Fig. 4b. 

Figure 4c shows how the spectral oscillations depend on the mean coronal electron 
temperature. As one would expect, the minimum in the relative amplitude, which 
corresponds to the point where the spectral oscillations go from their low-energy 
anticorrelation with the coronal density oscillations to their high-energy correlation, 
increases with coronal temperature. When Tinj is increased, both the phase lag and relative 
amplitude curves shift upward slightly in energy (Fig. 4d). 

4. Preliminary Application to Kilohertz QPO 

To illustrate how time-dependent Comptonization calculations can be utilised in data 
interpretation, we present an informal analysis of an 850 Hz QPO in 4U 1608-52, based 
on published observations of the energy-dependent QPO phase (Vaughan et al. 1997) and 
relative amplitude (Berger et al. 1996). With the assumption that oscillations in the soft 
photon source drive the QPO, Vaughan et al. estimated from the observed phase lags that 
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the size scale of the coronal gas is between a few kilometers and a few tens of kilometers. 

We model the observed spectral variability as a superposition of "basis variabilities" 
driven by oscillations (i) in the soft photon source, (ii) in the electron density of the 
Comptonizing corona, and (iii) in the temperature of the coronal electrons. This approach 
allows the possibility that the energy dependence of the observed oscillation phase 
might derive from the underlying phase relationships of the basis variabilities and their 
different energy dependences. For example, a photon injection oscillation might combine 
with small temperature oscillation that lags it by a quarter cycle. Since the size of the 
temperature-induced spectral oscillation grows with increasing photon energy, the resulting 
spectrum would have phase lags that increase steeply with energy, but are unrelated to the 
time it takes for photons to propagate from lower to higher energies. To safely interpret 
the phase lags, one should account consistently for the relative amplitude data and the 
time-averaged spectrum. It is also important that coronal emission dominates the spectral 
range under consideration. Contamination by (steady) noncoronal emission will reduce the 
relative amplitude at low photon energies below the model's predictions. 

Since time-averaged spectra for 4U 1608-52 during the QPO observations are not yet 
publicly available, the coronal parameters we derive must be taken as preliminary, subject 
to future revision. For example, we found that acceptable fits to the QPO observations 
could be produced with a wide range of coronal temperatures Tq. Here, for purposes 
of illustration, we assumed a temperature Tq = lOkeV, and an injection spectrum with 
^injo — 0.1 keV. We adjusted the remaining parameters by hand until an acceptable fit was 
obtained. Our "best fit" appears in Fig. 5. The coronal parameters for the fit are: I — 6 km 
and A'esco = 20 (or equivalently, t — A). The composition of the driving perturbation is 
"nei/neo = 0.32, Ni/Nq = — O.Sngi/neo, and Ti/Tq = — (0.37 + 0.08i)nei/neo- Although 
there is considerable latitude in acceptable admixtures of the injection and temperature 



- 17- 

perturbations, we found that to successfully match the relative amplitude data, the primary 
source of the QPO has to be an oscillation in the coronal electron density. Oscillations 
dominated by injection or coronal temperature oscillations produce unacceptable fits, even at 
the qualitative level. 

5. Conclusions 

Since the spectral variations driven by oscillations in the corona's supply of soft 
photons, heat, and material each bear their own distinctive characteristics, one can 
determine the nature of the oscillations driving QPO. Our sample application to a kilohertz 
QPO observed in 4U 1608-52 indicates that the main driver behind the QPO is an oscillation 
in coronal electron density. In agreement with earlier studies (Vaughan et al. 1997), we find 
that the size of the corona is comparable to estimates of the neutron star radius, supporting 
the notion that kilohertz QPO originate very close to the accreting neutron star. 
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Fig. 1. — When the Comptonization process is disturbed in different ways, distinctive 
patterns of spectral variabihty result. Shown are the spectral variations driven by 1% 
oscillations in (i) the soft photon injection rate (short dashed curves), (ii) the coronal electron 
temperature (unbroken curves), and (iii) the coronal electron density (long dashed curves). 
The driving oscillations are all at 40 Hz. Bold lines show the rms relative amplitudes of 
the spectral oscillations as functions of photon energy. An rms relative amplitude of 1% 
corresponds to a 1% oscillation in photon count rate about its average value. Light lines 
show the phase lag between the spectral oscillation and the physical oscillation driving it. 
Each spectral variation is produced in a "standard" corona with the following time-averaged 
parameters: A^esco = 100, Tlnj = 0.1 keV, Te = lOkeV, and Zgize = lO^cm. 

Fig. 2. — Spectral variations driven by 1% oscillations in the photon injection rate are 
illustrated in Figs. 2a-2d. For each variation, a heavy curve denotes the relative amphtude 
of the spectral oscillation, while the phase lag is depicted by a light curve. Each panel shows 
how the spectral variations depend on a particular parameter (the oscillation frequency, the 
average number of collisions suffered by an escaping photon, the coronal temperature, or 
the temperature of the injected photons), while the other parameters are kept at the values 
employed in Fig. 1. Fig. 2a compares the spectral responses at different driving oscillation 
frequencies. Variations induced by photon injection oscillations with frequencies of 40 Hz 
(unbroken curves), 100 Hz (short dashed curves), and 1 kHz (long dashed curves) are shown. 
Fig. 2b displays variations in the output spectra of coronae having Aesco — 10 (short dashed 
curves), 100 (unbroken curves), and 400 (long dashed curves). Fig. 2c shows variations in the 
output spectra of coronae with electron temperatures of 5 keV (short dashed curves), 10 keV 
(unbroken curves), and 15 keV (long dashed curves). Fig. 2d contains spectral variations 
with several different input photon temperatures: 0.1 keV (unbroken curves), 0.5 keV (short 
dashed curves), and 1 keV (long dashed curves). 
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Fig. 3. — Spectral variations driven by 1% oscillations in the coronal temperature. As in 
Fig. 2, each panel shows how the spectral oscillations depend on the underlying physical 
properties of the Comptonizing corona and its supply of soft photons. 



Fig. 4. — Spectral variations driven by 1% oscillations in the coronal electron density. Each 
panel shows how the spectral oscillations depend on the underlying physical properties of 
the Comptonizing corona and its supply of soft photons. The plotting conventions are the 
same as in Fig. 2. 



Fig. 5. — Sample fit to an 850 Hz QPO observed in 4U 1608-52. The data for the relative 
amphtudes are from Berger et al. 1996, and the phase data are from Vaughan et al. 1997. 
The relative amplitude data imply that the QPO is caused mainly by an oscillation in the 
coronal electron density, rather than by oscillations in the soft photon source or in the coronal 
temperature. Details of the other coronal parameters are given in the text. 
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